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Abstract 

Background: Vascular endothelial growth factor A (VEGFA) is a major regulator of both physiological and 
pathological angiogenesis. Associations between polymorphisms in VEGFA and complex disease have been 
inconsistent. The parent from whom the allele was inherited may account for these inconsistencies. This study 
examined the parent of origin effect on the expression of murine Vegfa. 

Methods: Two homozygous, inbred mouse strains A/J (AJ) and 129x1/SvJ (129) were crossed to produce reciprocal 
AJ129 and 129AJ offspring, respectively. RNA was extracted from cardiac tissue of 6 week old male (n = 8) and 
female (n = 8) parental, and male and female F1 offspring mice (AJ129 n = 8 and 129AJ n = 8). Vegfa and Hifla 
expression levels were measured by qPCR and compared between the F1 offspring from the reciprocal crosses. 

Results: We found significant differences in the expression of Vegfa in F1 offspring (AJ129 and 129AJ mice) of the 
reciprocal crosses between AJ and 129 mice. Offspring of male AJ mice had significantly higher expression of Vegfa 
than offspring of male 129 mice (p = 0.006). This difference in expression was not the result of preferential allele 
expression (allelic imbalance). Expression of Hifla, a transcriptional regulator of Vegfa expression, was also higher in 
F1 offspring of an AJ father (p = 0.004). 

Conclusion: Differences in Vegfa and Hifla gene expression are likely the result of an upstream angiogenic 
regulator gene that is influenced by the parent of origin. These results highlight the importance of including 
inheritance information, such as parent of origin, when undertaking allelic association studies. 
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Background 

Angiogenesis is the growth of capillary sized vessels. It 
plays a critical role in wound healing and embryological 
development, and in the pathogenesis of diseases such as 
cancer. Promoting angiogenesis can aid in the treatment 
of ischemia and myocardial infarction, while inhibition of 
angiogenesis has been effective in the treatment of cancer. 
Vascular endothelial growth factor A (VEGFA) is one of 
the most potent inducers of angiogenesis, and is a major 
regulator of both physiological and pathological angiogen- 
esis. It has been shown to stimulate blood vessel formation 
in animal models of ischemic cardiovascular disease [1], 
has been trialled as a potential therapeutic for myocardial 
ischemia [2], and therapies that inhibit VEGF are used in 
the treatment of cancer [3] . 
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A number of studies have investigated associations be- 
tween polymorphisms within, or nearby, the VEGFA 
gene and complex disease, including coronary artery dis- 
ease [4], atherosclerosis [5], myocardial infarction [6], 
and hypertension [7]. However, these sequence variants 
explain only a small fraction of the estimated heritability 
of the diseases studied, and the ability to replicate these 
findings can prove difficult. For example, inconsistent 
results have been shown when studying the association 
between the development of coronary artery lesions in 
individuals with Kawasaki disease with VEGFA polymor- 
phisms [8,9]. Population differences, study power, and 
choice of polymorphism may explain inconsistent re- 
sults. Additionally, the parent from whom the variant 
was inherited may also contribute to the inability to rep- 
licate results across studies. 

Parent of origin effects influence gene expression in 
offspring and result in the preferential expression of the 
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paternally or maternally inherited allele. It is caused by a 
variety of mechanisms including imprinting, sex, epigen- 
etic effects, and environmental influences during devel- 
opment [10]. In genetic studies, parent of origin effects 
are largely ignored; however, recent studies have found 
evidence of parent of origin effects in breast cancer [11], 
colorectal cancer [12], and type 2 diabetes [11]. These 
findings highlight the importance of investigating the ef- 
fects of parent of origin on phenotype and disease 
susceptibility. 

Due to the inconsistency in results when studying the 
association between VEGFA and disease, this study 
aimed to determine whether parent of origin effects in- 
fluence the expression of Vegfa. This study was carried 
out in cardiac tissue from Fl hybrids resulting from re- 
ciprocal crosses of two inbred mouse strains. 

Methods 

All experiments were carried out at the School of Medi- 
cine Animal Facility, University of Western Sydney and 
approved by the University of Western Sydney Animal 
Ethics Committee. Two homozygous, inbred mouse 
strains A/J (AJ) and 129xl/SvJ (129) were used. Crosses 
between AJ male and 129 female, and 129 male and AJ 
female produced reciprocal AJ129 and 129AJ offspring, 
respectively. Male (n = 8) and female (n = 8) mice from 
each of the parental strains were euthanised at six weeks 
of age by cervical dislocation. Male and female offspring 
were collected from each of the crosses (AJ129 n = 8 per 
sex and 129AJ n = 8 per sex). The hearts were excised, 
rinsed with PBS and flash frozen in LN2. 

Results 

Heterozygous Fl hybrids (AJ129 and 129AJ) resulting 
from reciprocal crosses of two inbred homozygous 
strains (AJ and 129) were used in this study. Direct se- 
quencing of the full Vegfa mRNA transcript in homozy- 
gous parental strains and heterozygous Fl offspring 
confirmed a single nucleotide polymorphism (SNP) at 



the 3'UTR of the gene (AJ = GG, 129 = AA, AJ129 and 
129AJ = GA). 

Vegfa expression 

This study compared the expression of Vegfa in AJ and 
129 parental mice, and in male Fl offspring from recip- 
rocal crosses (AJ129 and 129AJ), to determine whether 
parent of origin influences mRNA expression. Strain 
specific differences in Vegfa expression were seen be- 
tween AJ and 129 mice, with Vegfa expression 1.5 fold 
higher in male AJ mice compared to male 129 mice 
(Figure la). No strain specific differences in Vegfa ex- 
pression were observed between female AJ and female 
129 mice (Figure lb). In the male Fl offspring from the 
reciprocal crosses, Vegfa expression was 1.7 fold higher 
in AJ129 mice compared with 129AJ mice (Figure lc). 
To ascertain whether this difference in Vegfa expression 
tracked with Y chromosome transmission, the expression 
of Vegfa was measured in female parentals and offspring. 
There was no difference in Vegfa expression between 
parental AJ males versus females (p = 0.15), and female 
offspring of AJ males had 1.3 fold increased expression 
of Vegfa (p = 0.03) compared with females of the recip- 
rocal cross. 

Measurement of allelic expression 

To determine whether differences in Vegfa mRNA ex- 
pression between the Fl offspring from the reciprocal 
crosses were the result of unequal expression of the ma- 
ternally or paternally inherited allele, the expression of 
each allele was measured using quantitative ARMS PCR. 
The rsl6821436 SNP in the 3'UTR of Vegfa (G/A) 
was used to differentiate expression of alleles in the Fl 
heterozygous mice. In the Fl offspring, the expression of 
each allele relative to the homozygous parental strain 
was significantly lower, as expected. The expression 
of the G allele in the heterozygous AJ129 and 129AJ 
mice was significantly lower, 43% and 43%, than the 
homozygous AJ mice (Figure 2ai). Similarly, the ex- 
pression of the A allele in the heterozygous offspring 
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Figure 1 Relative expression of Vegfa. a Parental: AJ male compared to 129 male. b. Parental: AJ female com 
(offspring of AJ male) compared to 129AJ (offspring of 129 male). 
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Figure 2 l/egftj allelic expression, ai. Relative expression of the G allele in males, aii. Relative expression of the A allele in males, bi. Expression 
ratio of the G allele to the A allele in AJ129 mice. bii. Expression ratio of the G allele to the A allele in 129AJ mice. 



was significantly lower, 47% and 45%, than the homo- 
zygous 129 parental mice (Figure 2aii). No difference 
in the expression of the G allele relative to the A allele 
was seen within the Fl offspring (Figure 2bi and ii). 
Additionally, there was no difference in the expression 
ratio of the G allele to the A allele between the AJ129 
and 129AJ offspring. 

Hifla expression 

To determine whether differences in Vegfa expression 
were the result of trans-acting regulatory elements, the 
expression of Hifla, a transcription factor and known 
activator of Vegfa was measured. The expression of 
Hifla was measured in the AJ and 129 parental strains 
and in the Fl offspring of reciprocal crosses. The expres- 
sion of Hifla was 1.4 fold higher in male AJ compared 



to male 129 mice (Figure 3a). No difference in Hifla ex- 
pression was observed between the female parental 
strains (Figure 3b). In the male Fl mice, the expression 
of Hifla was 1.4 fold higher in AJ129 compared to 
129AJ mice (Figure 3c). 

Discussion 

This study examined the parent of origin effect on the 
expression of Vegfa, an important regulator of angiogen- 
esis. We found significant differences in the expression 
of Vegfa in Fl offspring (AJ129 and 129AJ mice) of the 
reciprocal crosses between AJ and 129 mice. Specifically, 
offspring of male AJ mice had significantly higher ex- 
pression of Vegfa than offspring of male 129 mice. This 
difference in expression was not the result of preferential 
allele expression (allelic imbalance), as no difference in 



2.00 
1.75 
1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 




2.00 
1.75 
1.50 
1.2S 
1.00 
0.75 
0.50 
0.25 
0.00 



II 



p = 0.004 



Figure 3 Relative expression of Hifla. a. Parental: AJ male compared to 129 male. b. Parental: AJ 
(offspring of AJ male) compared to 129AJ (offspring of 129 male). 
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the expression of the paternally or maternally inherited 
allele was observed in the Fl heterozygotes. Expression 
of Hifla, a transcriptional regulator of Vegfa expression, 
was also found to be significantly higher in Fl offspring 
of an AJ father. These results cannot be explained by 
Y chromosome transmission, as female offspring of AJ 
fathers had increased expression of Vegfa. In addition, 
the A} female parentals had no difference in expression 
of Vegfa when compared with the male AJ parentals. 
The DNA sequences of Fl heterozygotes are identical 
on autosomal chromosomes and differ only on their sex 
chromosomes. As Vegfa and Hifla are not located on 
sex chromosomes, this suggests that the observed differ- 
ences in Vegfa and Hifla gene expression are not a re- 
sult of DNA sequence variation and may be explained 
by a parent of origin effect. 

Angiogenesis is a strictly regulated process that results 
from a balance between pro and inhibitory angiogenic 
factors. VEGF is a potent pro-angiogenic factor and its 
expression is tightly controlled at a number of levels. 
Cz's-acting SNPs have been shown to alter Vegfa pro- 
moter activity [13-15], while SNPs in the 3'UTR have 
been reported to affect expression via the binding of mi- 
cro RNAs which affect mRNA stability [16,17]. The AJ 
and 129 parental strains differed at a single nucleotide in 
the 3'UTR of Vegf [18]. In this study, the paternally and 
maternally inherited alleles within the Fl offspring were 
expressed equally, indicating that ds-acting regulatory 
elements within, or nearby, the Vegfa gene were unlikely 
to account for the observed differences in Vegfa gene ex- 
pression. However, the presence of cis- acting regulatory 
elements in genes known to regulate Vegfa expression, 
which are influenced by parent of origin, may explain 
the difference in Vegfa expression. 

A number of tra«s-acting factors are known to regu- 
late the expression of Vegfa. Transcription factors that 
bind to response elements within the Vegfa promoter in- 
clude hypoxia-inducible factor (HIF1) [19], nuclear fac- 
tor (NF)-kB [20], activator protein (API) [21], and 
specificity protein (SP1) [22]. HIF1 is a transcription fac- 
tor that functions to restore oxygen homeostasis during 
hypoxia via the regulation of genes involved in angiogen- 
esis (Vegfa), glycolysis and vasodilation [23]. In this 
study, Hifla expression was significantly higher in the 
Fl offspring of AJ males which corresponded to in- 
creased Vegfa expression within these same animals. Im- 
balanced expression of alleles within the Hifla gene 
could not be assessed as the Fl mice were homozygous 
[18]. No studies have reported a role for epigenetic 
mechanisms in the regulation of Hifl expression how- 
ever, epigenetic regulation of von Hippel-Lindau (VHL), 
an ubiquitin ligase that targets HIFla for proteosomal 
destruction under normoxic conditions, has been 
reported. Hypermethylation of the VHL promoter in 



renal cell carcinoma has been correlated with increased 
HIFla expression [24], and deletion of miR-23b has been 
reported to reduce HIFla and VEGFA expression via its 
targeting of VHL [25]. Further work is needed to address 
whether epigenetic mechanisms are involved in the 
regulation of genes known to regulate Vegfa and Hifla 
expression in this mouse model. 

Transgenerational epigenetic inheritance is another 
mechanism that may explain the differences in expression 
patterns between the Fl offspring from the reciprocal 
crosses. It refers to genetic or environmental exposures in 
the parent, which affect patterns of gene expression in 
the offspring, without the offspring having inherited the 
allele or been exposed to the environmental agent [26]. 
Examples of transgenerational epigenetic inheritance in- 
clude DNA methylation, histone modifications, changes in 
chromatin structure and the action of small non-coding 
RNA all of which can result in epialleles. That is, alleles 
that are variably expressed in the absence of genetic het- 
erogeneity [27]. The epigenetic state is established during 
early embryogenesis and can be influenced by whether the 
allele is inherited from the mother or the father, as is the 
case with the agouti viable yellow allele (A vy ) in mice [28] . 
Parent of origin effects may be observed for a number of 
genes, if these genes belong to pathways that are ultim- 
ately regulated by epialleles. This may be the case in the 
present study, where a regulator of genes involved in 
angiogenesis is influenced by the parent of origin, and 
downstream this alters the expression of both Hifla and 
Vegfa. Further work is required to establish which genes 
are controlled by epigenetic mechanisms within the path- 
way controlling angiogenesis. 

Conclusions 

Our study demonstrates that parent of origin influences 
Vegfa expression in mouse cardiac tissue, and the gen- 
etic factors which determine gene expression are more 
complex than simple Mendelian inheritance. Addition- 
ally, given the common use of A/J and 129xl/SvJ strains 
to generate genetically engineered mice, studies using 
mixed genetic background mice, or studies using mice 
after backcrossing, need to consider parent of origin in 
addition to genetic variation. Finally, the identification of 
epigenetic regulation in genes involved in angiogenesis 
has important implications in how future genome-wide 
association studies need to incorporate technologies, 
which detect epigenetic and parent of origin informa- 
tion, to provide more informative associations. 

RNA extraction, quantification and cDNA synthesis 

Hearts were mechanically disrupted in the deep frozen 
state using the 'Mikro Dismembrator' automated frozen tis- 
sue disruptor (Sartorius, Germany), as previously described 
[29]. Total RNA was extracted from powderised tissue 
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using TRIzol (Sigma Aldrich, USA) as per manufacturer's 
protocol, followed by DNase I treatment using the 'DNA- 
free' kit (Applied Biosciences, USA) according to the man- 
ufacturer's protocol. RNA were quantified using the 
NanoPhotometer (Implen, Germany) and cDNA was 
synthesised from 500 ng of RNA using Bioscript™ Reverse 
Transcriptase and random hexamers (Bioline Pty Ltd, 
Australia) according to the manufacturer's protocol. Gen- 
omic DNA was extracted using the QIAamp DNA Mini 
Kit (Qiagen, Germany) as per the manufacturer's protocol. 

Quantitative PCR 

Quantitative PCR (qPCR) was used to measure mRNA 
expression of Vegfa and hypoxic inducible factor 1 alpha 
(Hifla). Mitochondrial ribosomal protein S28 {Mrps28S), 
18S ribosomal RNA (Rnl8S), and hypoxanthine guanine 
phosphoribosyl transferase {Hprtl) were used as normalising 
genes. Primer sequences are described in Table 1. Indi- 
vidual reactions (10 pi) contained 5x MyTaq PCR buf- 
fer (Bioline), 0.4 U of MyTaq (Bioline), lx SYBR Green 
I dye, forward and reverse primers (10 pM), and 4 ul of 
cDNA (or 4 pi water for no template controls). A 
pooled sample of cDNA was serially diluted 1 in 5 to 
generate a standard curve for each primer set, and the 
efficiency of the reaction was calculated. The PCR reac- 
tions were carried out in a MxPro3005P Real Time 
PCR System (Stratagene, Agilent Technologies, USA) 
under the conditions listed in Table 1, followed by a 
dissociation curve. All samples were run in triplicate. 

Direct sequencing 

Sanger sequencing was used to confirm the presence 
of a polymorphism in the 3'UTR of Vegfa (NCBI Refer- 
ence Sequence: NM_001025250.3). Total RNA from the 
parental strains and Fl offspring was amplified and 
the purified PCR products were sequenced (Australian 
Genome Research Facility, AUST). Sequence data was 
aligned and analysed using Sequencher v5.0.1 (Gene 
Codes Corp, MI). 



Measurement of allelic expression using amplification 
refractory mutation systems (ARMS) PCR 

Heterozygous Fl samples were evaluated for Vegfa allelic 
expression imbalance using quantitative ARMS PCR. 
The amplification primers (forward: CTTTCATCCCAT 
TGTCTACC(G/A); reverse: TGTTATTGGTGTCTTCA 
CTGGA) used in the assay generated a 188 bp amplicon. 
The allele specific primers were run in parallel reactions. 
Reactions were run under the following conditions: 95°C 
10 min; 95°C 30 sec, 62°C 30 sec, 72°C 30 sec for 30 - 
cycles; followed by a dissociation curve. Individual reac- 
tions (10 pi) contained 5x MyTaq PCR buffer (Bioline), 
0.4 U of MyTaq (Bioline), lx SYBR Green I dye, forward 
and reverse primers (10 pM), and 4 pi of cDNA (or 4 pi 
water for no template controls). All reactions were 
performed in triplicate. The ability of the primers to dis- 
tinguish between the two alleles with high specificity was 
assessed using cDNA from parental strains. Under these 
conditions neither primer exhibited cross-reactivity with 
the opposite allele. A standard curve was generated for 
each primer pair using a dilution series (200 ng to 0.20 ng, 
1 in 4) of heterozygous genomic DNA. This served as a 
reference for the 50:50 allelic ratio as we would expect a 
perfect 50:50 ratio of the two alleles in genomic DNA 
from a heterozygote with a diploid genome. We avoided 
using cDNA as a control as this would bias our results to- 
wards the allele expression ratio of the reference sample. 

Data analysis 

For qPCR, cycle threshold (Ct) values for each sample were 
calculated using the MxPRo QPCR software (Stratagene, 
Agilent Technologies, USA). Triplicate Ct values were aver- 
aged and the quantity (Q) of each sample was calculated 
using the delta-delta Ct method [30]. Q values from the nor- 
maliser genes were input into geNorm [31] and the geomet- 
ric means from these genes were used to generate a 
normalisation factor (NF). The Q values of the genes of 
interest were normalised by dividing by the NF value. Re- 
sults are expressed as mean ± SE fold changes in expression. 



Table 1 Primers and PCR conditions 


Primer 


bp 


Primer sequence 5' to 3' 


PCR conditions 


Hprtl F 


94 


TGACACTGGCAAAACAATGCA 


95°C 1 0' (95°C 40"; 60°C 30"; 72°C 30") X 40 


Hprtl R 


GGTCCmTCACCAGCAAGCT 


Dissociation curve 


Rn18SF 


178 


CTCTGGTTGCTCTGTGCAGT 


95°C 1 0' (95°C 40"; 62°C 30"; 72°C 30") X 40 


Rn18SR 


GGCTCCrrGTAGGGGTTCTC 


Dissociation curve 


Mrps28SF 
Mrps28SR 


223 


GAAATGCAAGCACGGAGAGT 
CCGATGACCAGTTTGTCCTT 


95°C 1 0' (95°C 40"; 62°C 30"; 72°C 30") X 40 
Dissociation curve 


VegfaF 


200 


GCTACTGCCGTCCGATTGAGAC 


95°C 1 0' (95°C 40"; 65°C 30"; 72°C 30") X 40 


VegfaR 


GTGCTGGCnTGGTGAGGTTTG 


Dissociation curve 


HiflaF 
HiflaR 


209 


GCTTCTGTTATGAGGCTCACC 
TCAAACTGAGTTAACCCCATGT 


95°C 1 0' (95°C 40"; 62°C 30"; 72°C 30") X 40 
Dissociation curve 
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Statistical significance between groups was determined 
using a general linear model. Statistical analysis was 
performed with PASW Statistics version 20. A P value 
of <0.01 was considered significant to account for 
multiple testing. 
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